clear;
t0=0;
tf=5;
y0=1;
f1=@(t,y)-(1.2+sin(10*t)).*y;
[t1,y1]=ode23(f1,[t0,tf],y0);
[t1,y1]

f2=@(t,y)cos(t)-y/(1+t^2);
[t2,y2]=ode45(f2,[t0,tf],1);
[t2,y2]

f3=@(t,y)[y(2);(2-2*t*y(2)-3*y(1))/(1+t*t)];
[t3,y3]=ode45(f3,[t0,tf],[0;1]);
[t3,y3]

f4=@(t,y)[y(2);y(3);cos(t)+5*cos(2*t)*y(3)/(t+1)^2-y(2)+y(1)/(3+sin(t))];
[t4,y4]=ode45(f4,[t0,tf],[1;0;2]);
[t4,y4]                                                                                                                  